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I. INTRODUCTION 

The extraordinary electric fields available from laser sys- 
tems make laser-driven charged particle acceleration an ex- 
citing possibility. Because of the maturity of RF acceler- 
ator technology, it is tempting to try to scale down tradi- 
tional microwave accelerators to optical wavelengths. How- 
ever, this proposal encounters several obstacles in the struc- 
ture design. First, because of the significant loss of metals at 
optical frequencies, we wish instead to use dielectric materials 
for a structure. Second, manufacturing a circular disk-loaded 
waveguide on such small scales poses a significant challenge. 
For these reasons, we must consider structures which differ 
significantly from those used in conventional accelerators. 

So far, macroscopic, far-field structures have been inves- 
tigated experimentally as a means to exploit high laser in- 
tensities for linear particle acceleration. In the Laser Elec- 
tron Acceleration Program (LEAP), a free-space mode was 
used to modulate an unbunched electron beam JT|. The ex- 
periment demonstrated the expected linear scaling of energy 
modulation with laser electric field as well as the expected 
polarization dependence. The next stage of the experimental 
program, to be performed at the E163 facility at SLAC, seeks 
to demonstrate net energy gain by first optically bunching the 
electron beam using an IFEL [2] [3). Beyond that, the desire 
is to demonstrate a scalable acceleration mechanism. For that 
purpose, a near-field, guided-mode structure would be more 
suitable. 

Photonic crystals [4| provide a means of guiding a speed- 
of-light optical mode in an all-dielectric structure. They have 
been investigated for some time for metallic RF accelerator 
structures because of their potential for eliminating a ma- 
jor source of beam break-up instability 16). In the opti- 
cal regime, a synchronous mode has been shown to exist in 
a photonic crystal fiber Q. Several years ago a study was 
conducted of two-dimensional planar structures (8). While 
the study was informative, such two-dimensional planar struc- 
tures are ultimately impractical because they only confine the 



accelerating mode in one transverse direction. In this paper we 
present a three-dimensional structure which overcomes this 
deficiency, confining the mode in both transverse dimensions. 
In Sec. Ill] we present the structure geometry and acceler- 



ating mode. Then in Sec. Ill we analyze the performance of 



the structure in terms of accelerating gradient and optical-to- 



beam efficiency. In Sec. IV we address single-particle beam 
dynamics in the structure, presenting a focusing-defocusing 
lattice and demonstrating stable beam propagation. 



II. STRUCTURE GEOMETRY AND ACCELERATING 
MODE 

The accelerator structure is based on the so-called "wood- 
pile" geometry, a well-established three-dimensional pho- 
tonic crystal lattice designed to provide a complete photonic 
bandgap in a structure with a straightforward fabrication pro- 
cess (9J. The lattice consists of layers of dielectric rods in 
vacuum, with the rods in each layer rotated 90° relative to the 
layer below and offset half a lattice period from the layer two 
below, as shown with the coordinate system in Fig. [T] We 
consider laser acceleration using a wavelength of 1550 nm, in 
the telecommunications band where many promising sources 
exist iflOll . At this wavelength silicon has a normalized per- 
mittivity of e r = e/eo = 12.1 ATI . We let a be the horizontal 
lattice constant, and take c - a V2~. As described in [9| the 
lattice has a face-centered cubic structure and exhibits an om- 
nidirectional bandgap. Based on those results, we take the rod 
width to be w = 0.28a. The ratio of the width of the bandgap 
to its center frequency is 18.7%. 

We form a waveguide by removing all dielectric material 
in a region which is rectangular in the transverse x and y di- 
mensions, and extends infinitely in the e-beam propagation 
direction z- In addition, there are two modifications: In or- 
der to avoid deflecting fields in the accelerating mode, we 
make the structure vertically symmetric by inverting the up- 
per half of the lattice so it is a vertical reflection of the lower 
half. Second, we extend the central bars into the waveguide 
by 0.124a on each side in order to suppress quadrupole fields. 
The consequences of this perturbation will be discussed fur- 
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ther in Sec. IV The geometry, with a defect waveguide intro- 
duced, is shown schematically in Fig.|2]and visually in Fig. [3] 
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FIG. 3: A symmetric waveguide. 



FIG. 1: A diagram of 8 layers (2 vertical periods) of the woodpile 
lattice. 
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FIG. 4: The accelerating field seen by a speed-of-light particle, av- 
eraged over a lattice period, normalized to the accelerating field on 
axis, shown with structure contours for a transverse slice at z = 0. 
The inner rectangle denotes the interface between free space and the 
absorbing boundary in the simulation. 
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FIG. 2: The geometry of a vertically symmetric waveguide structure. 



The inversion of half the lattice introduces a planar defect 
where the two halves meet, but this waveguide still supports 
a confined accelerating mode. Indeed, the mode is lossless 
to within the tolerance of the calculation, placing an upper 
bound on the loss of 0.48 dB/cm. Its fields are shown in 
Fig. |4] For this mode a/ A = 0.364, so using a 1550 nm 
source determines a = 565 nm. The individual rods are then 
158 nm wide by 200 nm tall. The mode was computed using 
the finite-difference time-domain method, including a uniaxial 
perfectly-matched layer to include losses IfTZl . 
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III. PERFORMANCE OF THE ACCELERATING MODE 

We can now explore the performance of this structure from 
the point of view of accelerating gradient and optical-to-beam 
efficiency. The gradient is limited by optical breakdown to the 
structure material. We can define a parameter, which is a prop- 
erty of the mode itself and independent of the material, which 
relates the accelerating gradient to the breakdown threshold 
of the dielectric. While the mechanism of optical breakdown 
is not fully understood for most materials, let us suppose that 
breakdown occurs when the electromagnetic energy density 
exceeds a certain threshold in the material. We then define 
the damage impedance of an accelerating mode by 

Z-^acc 
d = 7. . 

where £ acc is the accelerating gradient on axis and « max is 
the maximum electromagnetic energy density anywhere in the 
material. Then we can write the sustainable accelerating gra- 
dient as 

E. dcc = ^2Z d u th c. (1) 

For the mode in the woodpile waveguide, we have Zd = 
11.340. Optical breakdown studies in silicon have shown 
a breakdown threshold of M t h = 13.3 J/cm 3 at A = 1550 nm 
and lps FWHM pulse width [13|, resulting in an unloaded 
accelerating gradient of 301 MV/m. 

With designs for future high-energy colliders calling for 
beam power exceeding 10MW, optical-to-beam efficiency is 
a critical parameter for a photonic accelerator. Here we 
compute the efficiency, and the characteristics of the particle 
beam required to optimize the efficiency. Unlike conventional 
metallic disk-loaded waveguide structures, the waveguide we 
consider here has a group velocity a significant fraction of the 
speed of light; in this case v g = 0.269c. As the particle bunch 
and laser pulse traverse the structure, the bunch will gener- 
ate wakefields in the waveguide which destructively interfere 
with the incident laser pulse. Because the beam-driven wake- 
field has amplitude proportional to the bunch charge, wake- 
fields limit the practical amount of charge one can accelerate: 
If q is the bunch charge, then the energy gained by the bunch 
scales as q, but the energy loss due to wakefields scales as q 2 . 
Increasing the charge too much ultimately decreases the en- 
ergy gain, and thus the efficiency, of the structure. To improve 
the efficiency, we can embed the structure within an optical 
cavity in order to recycle the remaining laser pulse energy for 
acceleration of subsequent optical bunch trains. 

To compute the efficiency, we follow the treatment de- 
scribed in [14|. According to that description, the efficiency 
of the structure depends upon several parameters. First, the 
characteristic impedance of the mode, which describes the re- 
lationship between input laser power and accelerating gradient 
El, is 

E 2 r A 2 

z = _a £L _ =4840 
p 

where P is the laser power. Second, the group velocity of 
the mode affects the efficiency, as modes with group velocity 



closer to the speed of light couple better to a relativistic parti- 
cle beam. This is quantified by the loss factor, which is given 
bylEl 

= 1 cp g Z c 
41-/3 g A 2 ' 

where f3 g = v g /c. An optical bunch with charge q will radi- 
ate fields in the accelerating mode with decelerating gradient 
equal to kq. Finally, the Cerenkov impedance Z# parametrizes 
the energy loss due to wide-band Cerenkov radiation. Follow- 
ing ifTTl . we can estimate Z# from its value for a bulk dielec- 
tric with a circular hole of radius R, which is 

Z = z ° 

H 2n(R/A) 2 ' 

where Zq = 376.73 D. is the impedance of free space. For arbi- 
trary accelerating structures, we can let R be a length charac- 
terizing the radius of the vacuum waveguide, even if the guide 
is not circular. In our case the aperture is rectangular, so we 
define the parameter R by R = VA/2 = 0.451/1, where A is 
the aperture area. This yields Z# = 295 CI. 

Following 1 14 1, we assume the accelerator is inside an opti- 
cal cavity with a round-trip loss of 5%, and the laser pulses 
are coupled into the cavity through a beamsplitter with re- 
flectivity r. We consider a train of N optical bunches, each 
with the same charge, and assume that the duration of the train 
is much less than the slippage time At between the particles 
and the laser pulse inside the structure. We also assume that 
At = 3cr r , where cr T is the laser pulse duration. 

As described above, trying to accelerate too much charge 
can reduce the efficiency of an accelerator due to wakefield ef- 
fects. Therefore, for a given r and N, the structure has a max- 
imum efficiency ?7 max . For each N, we choose r to maximize 
^max- These optimum values are plotted in Fig. [5] From this 
we see that the efficiency can be made quite high. Even with 
just a single bunch, the efficiency reaches 38%, while for a 
train of 100 bunches, the efficiency is 76%. Once the optimum 
r is computed for each N, the total train charge q, and external 
laser pulse amplitude are chosen so that they together satisfy 
two conditions: First, that the peak unloaded gradient is equal 
to the maximum sustainable gradient in the structure, which 
we take to be 301 MV/m from the above discussion. Second, 
we require that the charge maximize the efficiency. We plot 
the optimum charge as a function of the number of bunches 
in Fig. [6] As expected from the results in IfHl . the optimum 
total charge is low, only 1.41 fC for a single bunch and 12.3 fC 
for 100 bunches. With the charge having been computed, we 
can then find the average initial gradient (the gradient experi- 
enced by the first optical bunch in the train) for each N. We 
can also compute the induced energy spread on the beam as 
the difference between the average accelerating gradient expe- 
rienced by the first and last bunches in the train, relative to the 
average initial gradient. These quantities are plotted in Fig. [7] 
From this plot we see that wakefields have a serious effect on 
the acceleration. For just a single bunch, the average initial 
gradient is reduced from 218 MV/m (which is reduced from 
301 MV/m due to the advancement of the electrons with re- 
spect to the Gaussian laser pulse envelope) to 197 MV/m. The 
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Number of optical bunches FIG. 7: The average initial gradient and induced energy spread as a 

function of the number of optical microbunches. 

FIG. 5: The optimum reflectivity and the corresponding maximum 
efficiency as a function of the number of optical microbunches in the 
bunch train. 

IV. PARTICLE BEAM DYNAMICS 
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FIG. 6: The optimum total charge of a bunch train as a function of 
the number of optical microbunches. 



minimum at around 7 optical bunches is due to two compet- 
ing effects. First, as the number of bunches increases from the 
single-bunch case, the total charge increases, generating more 
wakefields which are recycled within the cavity. But we see 
from Fig.|5]that as the number of bunches increases, the reflec- 
tivity decreases, so a smaller fraction of these wakefields are 
recycled. More concerning is the effect on the energy spread 
of the beam. Even with only two bunches, the spread in gra- 
dient is 6.2%, and with 5 bunches, the spread is 16.5%. 



As remarked above, we adjusted the geometry of the struc- 
ture because of beam dynamics considerations, namely to 
make the structure symmetric in order to suppress deflect- 
ing fields. Here we examine the particle beam dynamics in 
this structure in more detail. Two concerns immediately arise. 
First, the structure has an extremely small aperture, approxi- 
mately A in each transverse dimension, so confining the beam 
within the waveguide presents a serious challenge. Second, 
since the structure is not azimuthally symmetric, particles out 
of phase with the accelerating field will experience transverse 
forces. As it turns out, an investigation of the second problem 
will yield insight into the first: The optical fields in the struc- 
ture provide focusing forces which are quite strong compared 
to conventional quadrupole magnets and can be used to con- 
fine a beam within the waveguide. Therefore, we begin with 
a description of the transverse forces in the woodpile struc- 
ture. We then propose a focusing scheme, and compute the 
requirements on the phase space extent of the particle beam. 

In our analysis of the beam dynamics in this structure, we 
use the average of the fields over one longitudinal period of 
the photonic crystal waveguide. This is justified because the 
submicron period together with field intensities limited by 
damage threshold prohibits particles from acquiring relativis- 
tic momentum on the scale of a single period. To see this, one 
need only examine a normalized vector potential of the fields, 
defined by 

eEa 

A = — 2 ' 

mc L 

where E is the electric field magnitude. If A <K 1, then the 
variation of a particle's trajectory within a period is negligi- 
ble, so the time-averaged force on the particle is well approx- 
imated by the force due to the longitudinally averaged fields. 
In our case, E = 301 MV/m, so A = 3.3 x 10~ 4 . We therefore 
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proceed to define the longitudinally-averaged fields, as seen 
by a speed-of-light particle beam. Taking the fields to have 
e mt t j me dependence, we can define the averaged electric field 
amplitude by 



Average demodulated E normalized to E 



E(x,y) 



i r 

E(x,y,i 

a Jo 



z)e'^ z dz, 



where E is the complex electric field amplitude and k z = ko = 
a>/c is the Bloch wavenumber in the periodic structure. It fol- 
lows that the average amplitude of the transverse force on a 
particle is given by 



Thus E z essentially forms a "potential" for the transverse 
force. We note that the factor of i implies that particles ex- 
perience a transverse force to the extent that they are out of 
phase with the crest of the accelerating wave. 

Since the mode has speed-of-light phase velocity, E z sat- 
isfies V]_E Z = in the vacuum region. This implies that 
Vj. ■ Fj. = 0, so that focusing in both transverse directions si- 
multaneously is impossible. It also means that we can expand 
E z as a polynomial in the complex coordinate w = (x + iy)/A. 
Using the symmetries of the structure, we find we can write 
E, as 



E z = J]A 2m Re(w 2m ). 



(2) 



m=0 



We then have that 



dF* 
Ox 



= — A 2 . 

i=0,v=0 7tA 



Thus a particle ahead in phase by <p {kz — (ot = <f) will expe- 
rience focusing gradients K x and K v , in the x and y directions 
respectively, given by 



K x = -K y = ~Re — |< 
\7iAE 



eA 2 
~kAE ' 



(3) 



This is equivalent to a quadrupole magnet with focusing gradi- 
ent (A 2 /nAc) sin (p. The coefficients A 2m for m > 1 correspond 
to nonlinear forces. 

To examine the implications of these forces for our struc- 
ture, we extract the coefficients A 2m from the computed fields. 
We do so by fitting polynomials of a form similar to that in 
Eq.|2]for m = 0, . . . , 6 to the computed E x and E y fields; these 
are sufficient to compute the coefficients of the other compo- 
nents using the Maxwell equations. As noted in Sec. |ll] we 
have modified the structure geometry by extending the central 
bars into the waveguide in order to suppress the quadrupole 
fields. This insertion distance was determined by computing 
the fitted A 2 coefficient as a function of insertion distance, and 
finding the value for which A 2 = 0. This prevents particles 
from being driven out of the waveguide by the linear focusing 
force, but does not yet provide a mechanism for stable propa- 
gation within such a small waveguide. 
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FIG. 8: The accelerating field in the focusing structure modified to 
suppress the octupole moment. 



To address the question of stable propagation, we note from 
Eq.[3]that because of the optical scale of A, the focusing gra- 
dient can be quite large. This suggests using these strong fo- 
cusing forces to confine the particle beam by running a similar 
structure with the laser field 7r/2 out of phase with the parti- 
cle beam. We explore this possibility by again modifying our 
structure by inserting the central bars, this time to optimize 
the fields for focusing. For this purpose we wish to suppress 
the lowest order nonlinear field component, in this case the 
octupole field given by the A\ coefficient, for maximum sta- 
bility. We find that by inserting the central bars by 0.188a, 
we can suppress the octupole field, while the focusing gradi- 
ent for fields at damage threshold is equivalent to an 83 1 kT/m 
magnet. The accelerating field for this structure geometry is 
shown in Fig. [8] 

Having found suitable optical modes for both acceleration 
and focusing, we can now describe a system of focusing ele- 
ments designed to confine the particle beam. We consider a 
focusing-defocusing (F0D0) lattice 1 , in which each cell con- 
sists of a focusing segment of length Lf run n/2 behind in 
phase for focusing in the x direction, followed by an accelerat- 
ing segment of length L a , a focusing segment of length Lf run 
-k/2 behind in phase for defocusing in the x direction, and 
then another accelerating segment of length L a Ifl8l . We are 
free to choose the parameters Lf and L„. Let K be the peak fo- 
cusing gradient for the focusing segment, given by Eq.[3]with 
E the energy of the ideal particle. Assuming an initial ideal 
particle energy of 1 GeV, we find that if we run the focusing 
structure at damage threshold (with an on-axis accelerating 



1 For the remainder of this section, we refer to the F0D0 lattice as simply the 
"lattice," not to be confused with the photonic crystal lattice. 
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FIG. 9: Initial particle positions, with colors indicating whether, and 
when, each particle exited the waveguide aperture. The bounds of 
the plot correspond to the physical aperture of the structure. 



gradient of 133 MV/m), we have K = 2.49 x 10 5 rrT 2 . For our 
lattice we choose Lf = 0.2/ = 401 [xm, which gives a fo- 
cal length of / = 10.0 mm. We also choose the betatron phase 
advance per half-cell to be n/4; then L a - y2f = 14.2 mm. 
The maximum beta function value is then f5 F = 48.4 mm, the 
length of the unit cell is L c = 2(L a + Lf) = 29.1mm, and 
the betatron period is 4L C = 117 mm. As the particle energy 
increases, we consider these lengths to scale as 1 / *\[K °c VE- 

The lattice described in the previous paragraph is stable 
under linear betatron motion. However, nonlinear transverse 
forces in the accelerating and focusing modes may cause in- 
stabilities which limit the dynamic aperture of the lattice and 
constrain the allowable emittance of the particle beam. To de- 
termine the effect of these forces and compute the required 
emittance of a particle beam, we perform full simulations of 
particle trajectories. These simulations include the dynamics 
of the particles in all six phase-space coordinates, with field 
values taken from the polynomial fits. 

To compute the dynamic aperture of the lattice, we simulate 
particles with initial positions uniformly distributed through- 
out the waveguide aperture. The initial transverse momenta 
of the particles are taken to be 0, as are the initial deviations 
of the particle energies from the ideal case. We first simulate 
particles which are exactly on crest initially. We propagate 
the particles for 3 m through the lattice and record whether 
or not each particle collides with the waveguide edge. The 
initial particle positions, color-coded as to the particle's sur- 
vival, is shown in Fig. [9] In that figure, a red marker indicates 
that a particle initially at that position exited the waveguide 
aperture within 1 m of propagation. A green marker indicates 



that it remained within the waveguide through 1 m, but exited 
before 3 m, while a blue marker indicates that the particle re- 
mained in the waveguide through all 3 m of propagation. We 
see from the plot that most particles that are driven out of the 
waveguide exited within the first meter. The figure also indi- 
cates that the dynamic aperture extends to almost the entire 
physical aperture of the structure (the dynamic aperture being 
wider in x than y because the initial longitudinal particle posi- 
tions are at the high-/3 A point of the lattice). The initial energy 
of the particles is 1 GeV, and after 3 m the energy of the ideal 
particle is 1.87 GeV. The ellipse indicated in the figure is the 
ellipse of largest area which still contains only particles which 
remain in the waveguide. If we take this ellipse to be the 3cr 
boundary of the particle distribution, we can use the relation 
<x 2 = /3s to obtain the emittance requirement of the lattice. We 
find normalized emittances of 



9.2 x I0~ w m, 



„(«) 

y 



1.09 x 10 _9 m. 



(4) 



While these emittances are several orders of magnitude 
smaller than those produced by conventional RF injectors, be- 
cause of the small bunch charge as discussed in Sec. Ill the 



transverse brightness (charge per phase space area) required 
in the optical structure is similar to that of conventional accel- 
erators. 

We can also examine the beam dynamics for particles 
which are slightly out of phase with the crest of the accel- 
erating laser field. Running the same simulation as above, but 
with the particles given an initial phase offset, we find that the 
dynamic aperture is reduced as particles become further out 
of phase. This is shown in Fig. [10] For particles ahead in 
phase (4> < 0), we can see in the plots a pattern characteristic 
of a fourth-order resonance, indicating that the octupole field 
in the accelerating structure becomes significant as particles 
become off-crest. The difference in aperture shape between 
positive and negative phase shifts could be attributed to the 
sign difference of that octupole field. 

With simulations of the dynamic aperture of the lattice for 
both on- and off-crest particles, we can now compute the dy- 
namics of a realistic particle bunch, with variation in all six 
phase space coordinates. We take the initial normalized trans- 
verse emittances to be those computed above for on-crest par- 
ticles and given in Eq. [4] We also take the RMS phase spread 
and relative beam energy deviation to be cr$ = 10 mrad and 
cr a = 10~ 3 , respectively. We track the ensemble of particles 
for 3 m of propagation, as before, and record the phase space 
coordinates every 10 cm. We find that 261 of the 13228 parti- 
cles simulated, or 2.0%, exited the waveguide before complet- 
ing the 3 m propagation. From the recorded phase space coor- 
dinates we can track the normalized emittances as a function 
of propagated distance. The transverse emittances are shown 
in Fig. 
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We see that the emittances increase initially, but 
that this increase slows as the propagation continues. For the 
full 3 m, £•"' increases by 10.1%, while Sy grows by 7.3%. 
Thus the growth in invariant emittance is much smaller than 
the energy gain, so that the geometric emittance of the beam 
will still be adiabatically damped. This shows that particle 
bunches can be propagated stably in a photonic crystal accel- 
erator without an extraordinary enhancement to beam bright- 
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FIG. 10: Initial particle positions, with colors indicating whether, and when, each particle exited the waveguide aperture, for initial phases 
ranging from -90 mrad to 90 mrad. 



ness. 



V. CONCLUSION 

We have found numerically an accelerating mode in a three- 
dimensional dielectric photonic crystal structure. This mode 
propagates with speed-of-light phase velocity and is confined 
in both transverse directions with low loss, and it is there- 
fore suitable for acceleration of relativistic particles. Based 



on breakdown measurements for bulk dielectric, we estimate 
that the structure can sustain unloaded gradients in excess of 
300 MV/m. We computed that the structure can operate with 
high optical-to-beam efficiency using the mechanism of recy- 
cling laser energy in an optical cavity. However, the energy 
spread among optical bunches in a train remains a concern; 
some mechanism must be found to counter this effect or such 
an accelerator would need to operate with less than optimum 
charge. 

We also showed a novel method for beam confinement in a 
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FIG. 1 1 : Transverse emittances of the particle bunch as a function of 
propagated distance. 



small aperture which uses optical fields in a modified structure 
to focus the particle beam. We found that this mechanism only 
required brightness comparable to conventional accelerators, 
and resulted in minimal emittance growth even for significant 
propagation distance. 

Finally, while fabrication of these structures at optical 
lengthscales is likely to be challenging, much work on fabri- 
cating the woodpile lattice is underway in the photonics com- 
munity. Woodpile lattices with a small number of layers have 
been constructed using a layer-by-layer technique [ 19 1 and via 
wafer fusion |20|. Other experimental techniques, such as in- 
terferometric alignment GTl and stacking by micromanipula- 
tion l22ll hold promise for efficient, cost-effective fabrication. 
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